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SUMMARY 

Trajectories of water droplets about an ellipsoid of revolution 
with a fineness ratio of 10 (lO percent thick) in flight through a drop 
let field were computed with the aid of a differential analyzer. Analy 
ses of these trajectories indicate that the local concentration of' 
liquid water at various points about an ellipsoid varies considerably 
and under some conditions may be several times the free-stream concen- 
tration. Curves of the local concentration factor as a function of 
spatial position were obtained and are presented in terms of dimension- 
less parameters that describe flight and atmospheric conditions. The 
data indicate that the expected local concentration factors should be 
considered when choosing the location of devices that protrude into the 
stream from aircraft fuselages or missiles^ or when determining anti- 
icing heat requirements for the protection of these devices. 


INTRODUCTION 

The calculations presented herein were made in order to extend the 
study reported in reference 1 on the spatial distribution of cloud drop 
let concentration after the free-stream concentration and distribution 
are altered by an aerodynamic body moving through the cloud. In ref- 
erence 1 the alteration of the concentration of cloud droplets in the 
immediate vicinity of an ellipsoid of revolution with a fineness ratio 
of 5 (20 percent thick) is discussed; whereas, the data herein apply to 
an ellipsoid with a fineness ratio of 10 (lO percent thick). 

As is discussed in the reference cited, an aircraft moving through 
a cloud alters the concentration. of cloud droplets in the immediate 
vicinity of the aircraft. During flight in supercooled clouds, an un- 
even distribution of ice is often formed on rod-shaped objects that 
protrude from the fuselage. An illustration of this peculiar distribu- 
tion is shown in photographs presented in reference 1 and in figure 1 
of this report. The pitot mast shown in figure 1 was attached to the 
fuselage of a B-24 airplane. As indicated by this ice formation, there 
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is frequently a region of reduced (or zero)' droplet concentration next 
to the aircraft surface, followed by a narrow region of greatly in- 
creased droplet concentration farther out. Beyond the narrow region of 
high concentration, the droplet concentration gradually decreases toward 
the free -stream value with increase in distance from the surface of the 
aircraft. A sketch of a typical spatial variation is shown in figure 2. 

A knowledge of this spatial variation of local droplet concentra- 
tion about an aircraft or missile during flight throiigh clouds, drizzle, 
or rain is often important when choosing the location of devices that 
protrude into the stream or when determining the design requirements for 
icing protection on the devices. Examples of such devices are: (l) 
intake ducts and vents, (2) antenna masts, (3) ice detectors, and (4) 
instruments for measuring liquid-water content and droplet-size distri- 
bution. Similarly, at places where a body of revolution joins an air- 
foil (i.e., where there is a forebody ahead of the airfoil, such as a 
rocket pod or aircraft fuselage) , the local impingement of cloud drop- 
lets on the airfoil, and therefore the ice-protection design require- 
ments, will be altered by the effect of the body of revolution on the 
local droplet concentration. 

The nonuniform distribution of ice shown in figures 1 and 2 can 
be explained from a knowledgb of the trajectories of cloud droplets 
around the fuselage. Droplet trajectories about an ellipsoid of fine- 
ness ratio of 10 (lO percent; thick) in axisymmetric , inconq)ressible 
flow were calculated with t^e aid of a differential analyzer at the 
NACA Lewis laboratory. A pijolate ellipsoid is chosen as a good approx- 
imation for many fuselages and missile bodies, and has the additional 
advantage of a surrounding flow field that is known exactly for incom- 
pressible, nonviscous flowj A knowledge of the flow field is necessary 
for the computation of the/ droplet trajectories , which were calculated 
back to the midpoint of the ellipsoid and in space in the radial direc- 
tion as far out as a dista^ice equal to the minor axis. These trajec- 
tories were analyzed to determine the relation between the droplet con- 
centration at various points in space with the following variables: 
ellipsoid length and velocity, droplet size, flight altitude, and air 
temperature. The results of the analysis are stunmarized %n this report 
in terms of the dimensionless parameters ReQ and K, which include 

these variables. Although the calculations were made for incompres- 
sible flow, they should be applicable throiighout the subsonic region 
because of the small effect of compressibility on droplet trajectories 
(ref. 2) and because of the high flight critical Mach number of the 
ellipsoid. 


SYMBOLS 


The following symbols are used in this report: 
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ann\olar area perpendicular to major axis through which droplet 
fliix flows, sq ft 

local concentration factor, d(rQ)/d(r^), dimensionless 
droplet diameter, microns 

flux density of liquid water, lb/(hr)(sq ft) 

inertia parameter, 1.704X10“^^d^u/|iL, dimensionless (density 
of water, 1.94 slugs/cu ft, is included in constant) 

length of ellipsoid major axis, ft 

free-stream Reynolds number with respect to droplet diameter 
4.813X10“%p u/p, dimensionless 

free-stream Reynolds number based on volimie-median droplet 
diameter , dimensionless 

cylinder coordinates, ratio to major axis, dimensionless 
(fig. 3) 

starting ordinate at z = - <*> of droplet trajectory, ratio to 
major axis, dimensionless 

U free-stream velocity or flight speed, mph 

u ratio of local air velocity to free-stream velocity, 

dimensionless 

v' local droplet velocity, mph 

w' local liquid-water content, g/cu m 

p viscosity of air, slugs/(ft) (sec) 

density of air, sliags/cu ft 
Subscripts : 

av average of quantity over area A 

r radial conponent 

s ellipsoid surface 

z axial component 

0 free-stream conditions 


d 

F 

K 


Re 


0 


^®0,med 


r ,z 


rn 
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METHOD OF COMPUTING DROPLET TRAJECTORIES 

The eqmtions that describe the motion of cloud droplets about an 
ellipsoid are given in reference 3. A solution of the differential 
eqmtions of motion was obtained with the vise of a mechanical analog 
(described in ref. 4) based on the principle of a differential analyzer. 
The procedure for calculating the trajectories of cloud droplets with 
respect to the ellipsoid is the same as that described in references 1, 

3, and 5. As shown in figure 3, the ellipsoid orientation in the coordi- 
nate system used in references 1, 3, and 5 is retained herein. Since a 
flow field is axisymmetric around an ellipsoid of revolution oriented 
at 0° between its major axis and the direction of the free stream, the 
droplet trajectories in the elliptical section of all meridian planes 
are the same, and the spatial variation of concentration can be obtained 
from trajectory calculations in the z,r plane of figure 3. 

The equations describing the air-flow field around the body, re- 
quired for the solution of the equations of droplet motion, are given in 
reference 3. The values of the air-velocity components u^ and u^ as 

functions of r and z are given graphically in reference 5. 


RESULTS AND DISCUSSION OF TRAJECTORY COMPUTATIONS 
Dimensionless Parameters 

The equations of motion were solved for various values of the 
parameter i/k between 0.1 and 90. The Inertial parameter K is a 
measure of the droplet size, the flight speed and size of the body of 
revolution, and the viscosity of the air, in the form 

K = 1.704X10"^^ — (l) 

jUli 

The density of water, which is expressed as part of the conversion 
factor, is 1.94 slugs per cubic foot. For each value of the parameter 
i/K, a series of trajectories was computed for each of several values 
of free -stream Reynolds number ReQ. The free -stream Reynolds number 
is defined with respect to the droplet diameter as 

Reo = 4.813X10"® dp^U/p (2) 

In order to make these dimensionless parameters more physically sig- 
nificant in the following discussion, some typical combinations of K 
and Rsq are presented in table I for various lengths and velocities 

of the ellipsoid, droplet sizes, and flight pressure altitudes and 
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temperatures. A proced\are for rapid calculation of i/k and Roq from 
practical flight conditions is given in appendix B of reference 6. 


Average Mass Flux Density of Water in Droplet Form 

The average mass flux of water in droplet form per unit area through 
an annular area of space (of width - r^,' fig. 3) perpendicular to 

the major axis of the ellipsoid is obtained from the law of conservation 
of matter. Assume that the liquid water in droplet form is moving be- 
tween two s\irfaces formed by rotating two neighboring trajectories in 
the r^z plane about the axis of the ellipsoid^ as shown in figure 3. 

Then, 

“ '^av '^z,av ^ 

and 

^av _ ^av^z ,av _ ^0 ^ 0,2 ~ ’^O,! ^(^0^ 

Fo “ wqU ~ a ~ r2 - r2 " A(r^ 

^ X 

( 

The subscript 0 refers to conditions at large distances ahead of 
the ellipsoid (free-stream conditions), and the subscript av refers 
to the average of a quantity over the annular area A. The annxalar area 
A is measured in a plane perpendicular to the major axis of the ellip- 
soid (fig. 3). From equation (4), the average flux density thro ugh an 
annular area A or a sector of the annular area A can be written as 
follows ; 


( 3 ) 

( 4 ) 


F 


av 


0.33wqU 


A(r^)’ 


lb/(hr)(sq ft) 


( 5 ) 


The constant 0.33 is a conversion factor for the units used. 

Curves of r^ as a function of r^ are presented in figure 4 for 

various values of Rsq and i/k at three constant z-positions, the 
midpoint (z = O), a point l/4 major axis length from the noze (z = -0.25), 
and the tip of the nose (z = -0.5). These curves were obtained from 
calciilated droplet trajectories. Equation (5) and fig\ire 4 may be used, 
for example, to determine the amount of water in droplet form passing 
in a spatial region of cross-sectional annular area A located outside 
the ellipsoid a distance r from the major axis. 
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Local Liquid-Water Flux Density 

The local mass fl-ux of water in droplet form per unit area perpen- 
dicular to the major axis at a point in the vicinity of the moving 
ellipsoid can he obtained from equation (s) by letting A(r^) approach 
zero as a limit. Then 

F = O.SSwqUC, lb/(hr)(sq ft) (S) 


where 


C = d(r2)/d(r2) 

The local concentration factor C at a point in space is obtained from 
the slope of the curves of r^ as a function of r^ (fig. 4) at the 

point of interest for the ReQ and K combination being considered. 

The concentration factor C as a function of r for selected 
values of Roq and K at z = 0, -0.25, and -0.5 is given in figure 

5. As was stated for the 20-percent -thick ellipsoid in reference 1, 
the concentration factor approaches 1 (free-stream value) at large 
values of r for all z-positions. As in the case of the 20-percent- 
thick ellipsoid, the concentration factor is less than the free-stream 
value (C < l) at the nose of the ellipsoid (z = -0.50) near the major 
axis (small values of r, fig. 5(m)). At z = -0.50 no difference was 
found, within the accxoracy of the computation, when Rcq was changed 

and the value of l/K was retained. 

The peak values of the concentration factor given in figure 5 are 
plotted in figure 6 as a function of i/k for various values of Rsq. 

The peak values are given only for the midpoint (z = O) and the quarter 
point (z = -0.25) on the ellipsoid, because peak values do not exist at 
the nose (z = -0.50). The peak values of C were not well established 
by the trajectory computations for all values of Roq and i/k. For 
some cases, such as those given in figure 5(b), the peak values might 
be 50 percent higher than shown. The values shown in figures 5 and 6 
are weighted for consistency with other data available on this subject, 
including the data of reference 1. 

The r-position of the peak value of C at z = 0 and z = -0.25 
is plotted as a function of i/k for constant Rcq in figure 7 . With 

the use of figures 6 and 7, the peak value of C and its r position 
can be determined for the value of Roq and i/k of interest. 
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The local liquid-water content in grams per cubic meter can be ob- 
tained at any point in space by dividing the droplet flux by the z- 
component of the local droplet velocity at that point: 

WqU 

w = — C, g/cu m (7) 


The value of the z-component of the droplet velocity v' was noted 

z 

during the calculations of the droplet trajectories. A value of 

- 1 C±1 percent) was observed between z = -0.40 and z = 0. A 

value of v^/u == 0.98 (±2 percent) was observed between z = -0.50 and 
z = -0.40. 


Shadow Zone 

The region of zero concentration adjacent to the surface of the 
ellipsoid, which is illustrated in figure 2 and is evident in figure 5, 
is called the shadow zone. This region is protected from droplet pene- 
tration by the air -flow characteristics ahead of and in the vicinity 
of the forward positions of the ellipsoid. The thickness (r - rg) of 

the shadow zone at each z-position of the body for various Rsq and K 

values is given in figure 8. Generally, the thickness of the shadow 
zone increases as z approaches 0, that is, from nose to midpoint of 
the ellipsoid. The shadow zone is 0 at the nose for all values of 
Reo and 1/K shown in figvire 8 and becomes of finite thickness at a 
z-position that depends on Roq and i/k. 


Comparison of 10- and 20-Percent -Thick Ellipsoids 

Because a conplete conparison of the 10- and 20-percent -thick 
ellipsoids that covers a wide range of flight and atmospheric condi- 
tions is beyond the scope of this report, the following limited com- 
parison is made for a set of conditions that occur rather frequently. 
However, a danger must be noted arising from the choice of only a few 
sets of conditions for conparison: the relative inportance of the 

different factors being conpared may change with other conditions not 
discussed. 

As might be expected, the shapes of the cxirves for peak values of 
concentration factor C and for shadow zone (figs. 6 and 8, respec- 
tively) resemble the curves describing the same factors in reference 1 
for the 20-percent-thick ellipsoid. The principal difference between 
the curves for the two ellipsoids is that the effect of concentration 
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and the extent of shadow are not as great for the 10-percent -thick 
ellipsoid as for the 2 0-percent -thick body. The effect of a change in 
body thickness on peak value of concentration factor C is shown in 
figure 9 for an ellipsoid 10 feet long moving at 300 miles per hour at 
15^000 feet altitude through droplets 20 microns in diameter. The ef- 
fect of change in body thickness on the shadow zone under the same con- 
ditions is shown in figure 10. The thin body deflects the droplet 
trajectories from straight lines less than the thick body. For a body 
of zero thickness^ the concentration factor would equal the free-stream 
value everywhere in space and the shadow zone is eliminated. 

The variation of local concentration factor with spatial position 
and Rcq and i/k is discussed in reference 1. The concentration 

factor is very sensitive to small changes in 1/K and ReQ, as can be 

seen in figure 5. The variation of concentration of liquid-water con- 
tent with z- and r-position and i/k and Rcq indicates that care 

should also be exercised with the 10-percent -thick body in locating 
instrument sensing elements and small inlets or vents that protrude 
from the surface and are sensitive to impinging water or ice formation. 

The values presented in figures 4 to 8 apply directly only to 
flights in clouds con5>osed of droplets that are all uniform in size. 

The effect of droplet-size distribution on concentration factor is dis- 
cussed in reference 1 with the use of several illustrative examples. 

The manner in which the concentration factor applicable to a 10-percent- 
thick ellipsoid is altered by droplet-size distribution is given here 
for the same conditions used in one of the illustrative exanqples in 
reference 1. 

With the use of the same flight and atmospheric conditions, a com- 
peirison with the 20-percent -thick ellipsoid is possible, as shown in 
figure 11. The abscissa is the distance from the ellipsoid surface, in 
inches, at z = 0. For both the 10- and 20-percent -thick ellipsoids, 
the major axis length is 20 feet, the flight speed is 300 miles per 
hour, the altitude is 15,000 feet, and the volume -median droplet diam- 
eter is 20 microns. For both bodies the peak value of the concentra- 
tion factor considerably decreases as the distribution broadens from 
uniform droplets to a "C" distribution (defined in ref. l). Also, the 
region that was a former shadow zone with zero concentration for the 
cloud condition with imiform droplets becomes an area with a small con- 
centration factor. A comparison of concentration-factor curves calcu- 
lated for droplet-size distributions estimated to be present in clouds, 
such as shown in figure 11, shows good qualitative agreement with the 
shape of actual ice formations obtained on rod-shaped objects in flight 
throtigh the clouds (fig. l) . 
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CONCLUDING REMRKS 

The significance of the spatial variation of liquid-water content 
and droplet size is discussed thoroughly in reference 1. The effect of 
the spatial variation on measiaring instruments and on the rotating 
multicylinder method is also evaluated in reference 1. The importance 
of avoiding measurements of liquid water and droplet size in the region 
of variation of droplet concentration in the vicinity of an aircraft 
fuselage is emphasized, because the correction of the indicated water 
content and droplet-size distributions to free-stream conditions would 
be very involved. 

The data of this report are applicable in a quantitative manner 
only to ellipsoids of revolution with a fineness ratio of 10. As was 
stated in reference 1 for a thicker body (fineness ratio of 5), these 
data also apply approximately in the vicinity of the nose section of a 
body of revolution that can be physically matched with the nose section 
of a 10-percent -thick ellipsoid of a given length. The data are par- 
tic\ilarly applicable if the afterbody does not appreciably alter the 
sir-flow field ahead of the ellipsoid section from that attendant to an 
ellipsoid. Because many bodies of revolution of interest are different 
in shape from an ellipsoid, the data of this report and reference 1 are 
primarily useful in pointing out in a q-ualitative manner the type of 
variation in liquid-water concentration that might be expected in the 
vicinity of a body of revolution in flight throiogh a droplet field. 


Lewis Flight Prop\ilsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, February 1, 1955 
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Figure 1. - Ice formation on pitot mast on B-24 aircraft after 
flight through supercooled clouds. 







Figure 3. - Coordinate system for droplet trajectory calculations. 
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Figure 4. - Continued. Square of starting ordinate as function of square of droplet ordinate at constant z-posltlon. 
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Figure 4. - Continued. Square of starting ordinate as function of square of droplet ordinate at constant z-posltlon 
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Figure 4. • Continued. Square of starting ordinate as function of square of droplet ordinate at constant z-posltlon. 
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Figure 4. - Concluded. Square of starting ordinate as function of square of droplet ordinate 
at constant z-posltlon. 
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Figure 5. - Continued. Variation of local concentration factor with radial distance r at constant axial 
position z. 
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Figure 5. - Continued. Variation of local concentration factor with radial distance r at constant axial 
position z. 







Figure 5. - Continued. Variation of local concentration factor with radial distance r at constant axial position 
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Figure 5. - Continued. Variation of local concentration factor with radial distance r at constant axial position 
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Figure 5. - Continued. Variation of local concentration factor with radial distance r at constant axial position 
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Figure 6. - Concluded. Variation of peak value of concentration factor with l/K for constant free- 
stream Reynolds number. 





Figure 7. - Variation of peak position of concentration factor with l/K for constant free-stream 
Reynolds number. 
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Figure 8. - Continued. Thickness (r - r^) of shadow zone as function of z-posltlon 
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(d) Free-stream Reynolds number, 1024. 

Figure 8. - Continued. Thickness (r - rg) of shadow zone as function of z-posltlon 






(e) Pree-3tream Reynolds number, 4096. 

Figure 8. - Concluded. Thickness (r - Tq) of shadow zone as function of z-posltlon 
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Figure 9. - Variation of peak value of concentration factor 
vlth ellipsoid thickness. Major axis, 10 feet; speed, 300 
miles per hour; droplet diameter, 20 microns; altitude, 
15,000 feet; Rbq, 1287; 1/K, 16.6; z = 0. 
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Figure 10. - Variation of shadow-zone thickness with ellipsoid 
thickness. Major axis, 10 feet; speed, 300 miles per hour; 
droplet diameter, 20 microns; altitude, 15,000 feet; RSq, 

1287; 1/K, 16.6; z = 0. 
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Figure 11, - Effect of droplet-size distribution on concentration factor for ellipsoid. Major 
axis length, 20 feet; velocity, 300 miles per hour; altitude, 15,000 feet; air temperature, 

1 F; volume-median droplet diameter, 20 microns 128.7; 33.1. 




